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Abstract 

Effect  of  interface  layers  on  wave  propagation  in 
a  composite  medium  has  been  investigated  in  this 
paper.  Two  systems  are  considered:  scattering  of 
waves  by  an  inclusion  with  a  thin  interface  layer 
and  guided  waves  in  a  laminated  plate  with  thin 
interface  layers.  In  the  first  case  the  presence  of 
interface  layers  is  found  to  modify  the  scattering 
cross-section  significantly.  In  the  second  case  it  is 
found  that  guided  waves  in  a  laminated  plate  are 
slowed  down  considerably  in  the  presence  of  soft 
interface  layers.  The  results  of  this  study  can  be 
used  to  characterize  interface  layers  by  ultrasonic 
techniques. 


1.  Introduction 

In  composites  (particularly  metal  matrix  com¬ 
posites  reinforced  by  fibers  or  particles),  it  is 
often  the  case  that  there  is  an  interface  layer  sur¬ 
rounding  the  fibers  or  particles  induced  by  pro¬ 
cessing  conditions.  This  layer  has  properties 
different  from  those  of  the  matrix  and  the  rein¬ 
forcing  phase.  The  strength  and  fracture  behavior 
of  the  composite  is  significantly  influenced  by  this 
interface  layer.  It  is,  however,  difficult  to  charac¬ 
terize  this  layer  non-destructively.  Ultrasonic 
techniques  provide  a  means  to  characterize 
mechanical  properties  of  the  interface  regions.  In 
this  paper  we  present  results  showing  the  effect  of 
interface  layers  on  scattering  cross  sections  of 
spheroidal  inclusions.  The  null-field  or  T  matrix 
method  [1,  2]  has  been  used  to  study  the  scatter¬ 
ing  problem.  The  interface  layer  is  modelled  as  a 
thin  shell,  the  equations  of  motion  of  which  enter 
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the  boundary  conditions  on  the  surface  of  the 
inclusion.  In  an  earlier  study,  Mai  and  Bose  [3] 
considered  the  effect  of  thin  viscous  layers  on 
long-wavelength  scattering  by  spherical  inclu¬ 
sions.  More  recently  Datta  et  al.  [4,  5)  have 
studied  scattering  by  spherical  inclusions  sur¬ 
rounded  by  thin  interface  elastic  layers.  In  both 
these  studies  it  is  assumed  that  the  tractions  are 
continuous  across  the  layer,  whereas  displace¬ 
ments  satisfy  jump  conditions  that  are  linear  in 
the  thickness  of  the  layer.  The  resulting 
approximate  boundary  conditions  are  based  on 
the  assumption  that  inertial  and  curvature  effects 
are  negligible.  The  present  study  shows  that  both 
these  have  a  strong  effect  on  the  scattering 
behavior. 

As  another  example  of  the  influence  of  inter¬ 
face  layers  on  wave  propagation  in  composites  we 
have  examined  guided  wave  propagation  in  a 
composite  laminate  with  interface  bond  layers.  A 
stiffness  method  [6,  7)  is  used  to  study  the  disper¬ 
sion  of  guided  waves  with  or  without  bond  layers. 
It  is  found  that  the  dispersion  depends  strongly 
on  bond  layer  properties. 

2.  Scattering  by  an  inclusion 

Consider  an  inclusion  with  a  thin  interface 
layer  of  thickness  h,  embedded  in  an  elastic 
matrix.  Let  S  be  the  surface  of  the  inclusion  with 
its  outward  pointing  unit  normal  n.  The  layer  is 
assumed  to  be  sufficiently  thin  that  all  boundary 
conditions  can  be  applied  on  S.  The  density  and 
Lame  constants  of  the  inclusion  are  taken  to  be 
p,,  2,  and  p,  respectively.  The  corresponding 
quantities  for  the  layer  and  the  matrix  are  p,„  A,„ 
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/u{)  and  p,X,p  respectively.  A  time  dependence  of 
the  form  exp(  -  iwr)  is  assumed  throughout. 

In  the  null-field  approach,  all  pertinent  fields 
are  expanded  in  spherical  vector  waves.  For 
brevity  we  present  here  only  an  outline  of  the 
method.  The  details  can  be  found  elsewhere  [1,  8, 
9j.  In  the  matrix  the  incident  regular  waves  are 
denoted  by  Re  ipn  and  the  outgoing  scattered 
waves  by  ip„,  where  n  is  the  quadruple  index.  The 
corresponding  transmitted  waves  in  the  inclusion 
are  Re  V',,'.  The  incident  and  scattered  waves  are 
then  expressed  as 

u"  =  Z  a„  Re  xpn  ( 1 ) 

n 

(2) 


The  problem  is  now  to  express  u „  and  / ^  in 
eqns.  (4)  and  (5)  in  terms  of  «_  and  I  and  to  this 
end  the  properties  of  the  layer  must  be  taken  into 
account.  As  h  is  assumed  small  (compared  with 
both  the  wavelengths  and  the  radii  of  curvature 
of  S)  a  consistent  approach  is  to  start  from  the 
three-dimensional  equations  for  the  layer  and 
approximate  to  order  h  wherever  possible  to 
obtain  a  theory  that  is  exact  to  order  h.  For  a 
rotationally  symmetric  inclusion  this  leads  to 


/  du\ 

\9«/.an 

_  (  f  -  )lan 

Ad 

-«x(V,x«_) 

(8) 

du 

n-  = 

1 

(n  t.  -x,ys-u. ) 

(9) 

3  n 

2d  +  2  ju„ 

Here  a„  is  assumed  known  and  f,  is  to  be  deter¬ 
mined,  i.e.  the  problem  is  to  determine  the  trans¬ 
ition  matrix  elements  Tnn  that  give  the  relation¬ 
ship  between  an  and  f, : 

f,=  Y.  Tnnan.  (3) 

ri 

Following  the  usual  procedure  in  the  null-field 
approach  the  outer  integral  representation  gives 
relations  between  n„,  /„  and  the  surface  displace¬ 
ment  and  traction  U  : 

««=  -  —  f  *+‘Vn)dS  (4) 

f*  s 

/»“—  f  («+-f(Re  ?„)-*♦ -Re  VOd 5  (5) 

f*  s 

Here  k s  is  the  shear  wavenumber  in  the  matrix. 

The  index  +  denotes  the  limit  from  the  out¬ 
side  and  t  is  the  traction  operator.  The  surface 
field  on  the  inclusion  is  expanded  in  the  regular 
spherical  waves  of  the  inclusion 

«-=!«„  Re  (6) 

and  the  inner  integral  representation  then  gives 
an  expansion  of  the  surface  traction  on  the  inclu¬ 
sion  as 

f-  =  Za/(  ReVC)  (7) 

n 

where  /'  is  the  traction  operator  of  the  inclusion. 


and 

(10) 

on 

The  expression  for  B  is  lengthy  and  can  be 
found  in  ref.  10.  It  contains  terms  representing 
the  inertial  and  curvature  effects.  It  should  be 
noted  that  a  most  common  assumption  that  is 
made  dealing  with  thin  layers  [5]  is  that  the 
normal  traction  is  continuous  across  the  layer. 
This  means  that  B  is  set  equal  to  zero.  Also, 
curvature  effects  represented  by  the  last  two 
terms  in  (8)  and  (9)  are  dropped. 

Using  equations  (4)— (6 )  it  is  now  possible  to 
derive  the  transition  matrix  elements  T„„.  This 
is  rather  complicated  and  will  be  omitted  here. 
Because  of  the  complexity  of  using  the  exact  O(h) 
approximation  to  Tnn •,  it  seems  worthwhile  to 
look  for  some  simplifications.  As  mentioned 
before,  a  commonly  made  simplifying  assumption 
would  be  to  neglect  all  the  inertial  and  curvature 
effects  in  ( 8)— ( 10).  Another  simplification  would 
be  to  use  a  membrane-shell-type  approximation 
for  the  thin  layer.  The  latter  simplifies  the  expres¬ 
sion  for  B  to 

B'(u,  t)=  - p{)aj2u+\so\u)  (11) 

This  expression  can  be  simplified  further  by- 
keeping  only  the  first  term  on  the  right-hand  side 
which  represents  the  inertial  effect.  This  gives 

B"=-pt](o:u  (12) 

In  ref.  9,  results  were  presented  for  the  scattering 
by  a  spherical  inclusion  with  a  thin  surrounding 
layer  using  the  exact  solution  to  the  problem  and 
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the  approximations  (,11)  or  (12).  It  was  found  that 
use  of  (12)  gave  very  good  results  at  low  frequen¬ 
cies  and  the  predictions  using  (11)  and  (12) 
started  to  diverge  at  high  frequencies.  In  Section 
4  we  present  some  results  for  the  scattering  of  P 
waves  by  an  SiC  inclusion  with  interface  layer  in 
an  aluminum  matrix,  showing  the  comparison 
between  the  predictions  of  the  full  O (h)  approxi¬ 
mation  and  the  three  simplifying  approximations 
mentioned  above. 

3.  Rayleigh -Lamb  waves  in  a  laminated  plate 

Consider  a  laminated  plate  where  each  lamina 
is  made  up  of  aligned  continuous  fiber  reinforced 
material  and  there  is  a  thin  bonding  material  layer 
between  adjacent  laminae.  In  this  paper  we  con¬ 
sider  a  cross-ply  plate  (0o/907--*/9070o)  and 
assume  that  the  waves  are  propagating  either  in 
the  0°  or  in  the  90°  direction. 

On  the  assumption  that  the  wavelength  is  much 
larger  than  the  fiber  diameter  and  the  spacing 
between  the  fibers,  each  lamina  can  be  modelled 
as  a  transversely  isotropic  medium  with  the 
symmetry  axis  aligned  with  the  fiber  direction.  In 
that  case  it  is  possible  to  derive  an  exact  disper¬ 
sion  equation  for  propagation  of  harmonic  waves 
in  the  plane  of  the  plate.  However,  for  a  larger 
number  of  laminations  this  equation  is  very  com¬ 
plicated  and  requires  extreme  care  in  locating  the 
roots.  For  this  reason,  we  have  presented  before 
[6.  7,  1  1)  a  stiffness  method  that  applies  to  an 
arbitrary  number  of  layers. 

In  this  method  each  lamina  is  divided  into 
several  sublayers.  In  each  sublayer  the  displace¬ 
ment  is  assumed  to  have  the  form 

«(.v.  y.  z,  t)=  U(z)  exp)i(  kx  +  /y  -  (Dt)\  ( 13) 

where  we  may  choose  the  x  axis  parallel  to  the  0° 
fibers  and  the  z  axis  in  the  thickness  direction  of 
the  lamina.  Choosing  the  origin  in  the  midplane 
of  the  sublayer,  L'\z)  is  approximated  by  poly¬ 
nomial  interpolation  functions  with  coefficients 


as  the  displacements  and  tractions  at  the  top  and 
bottom  of  the  sublayer.  This  assures  continuity  of 
these  quantities  at  the  interfaces  between  the  sub¬ 
layers  and  the  laminae.  It  has  been  shown  [6,  7, 

1  1  j  that  this  method  (method  I)  gives  very  accur¬ 
ate  results  even  at  high  frequencies.  However,  this 
method  leads  to  a  more  complicated  eigenvalue 
(algebraic)  problem  than  if  only  displacement 
continuity  at  the  interfaces  is  enforced.  In  the 
latter  case,  one  loses  the  high  accuracy  that  is 
obtained  using  method  I.  A  compromise  between 
the  two  is  that  used  in  refs.  12  and  13.  called 
method  II,  w'here  the  interpolation  polynomials 
involve  displacements  at  the  top,  middle  and 
bottom  surfaces.  It  has  been  shown  [6]  that  this 
last  method  gives  good  results  except  at  high  fre¬ 
quencies.  Because  method  II  is  somewhat  simpler 
to  use  when  a  large  number  of  laminae  has  to  be 
considered  we  have  used  that  in  this  paper.  Since 
both  method  I  and  method  II  have  been  well 
documented  before,  we  shall  omit  the  details 
here.  In  Section  4  we  present  numerical  results 
showing  the  effect  of  interfaces  on  the  dispersion 
behavior. 

4.  Numerical  results  and  discussion 

4. 1.  Effect  of  interface  layer  on  scattering 

In  this  section  we  present  some  results  for  the 
scattering  cross-section  of  spheroids.  We  con¬ 
sider  incident  longitudinal  waves.  As  an  example 
we  take  a  SiC  particle  embedded  in  an  aluminum 
matrix.  The  material  properties  of  the  matrix  are 

A  +  2//  =  1 10.5  GPa 

p  =  26.7  GPa 

p  =  2.706  kg  nT  ’ 

Those  of  the  particle  are 

A,  +  2^,  =474.2  GPa 

px  =  188.1  GPa 

pt  =  3.181  kg  nr  1 


TABLE  1  Comparison  of  analytical  and  numerical  results  for  the  total  scattering  cross-section  for  a  plane  P  wave  incident  on 
an  SiC  sphere  in  aluminum 


ksa 

Analytical 

It  =  t).  la 

f  ull  order 
h 

Membrane 

Inertial 

Only  jump 
in  u 

Order  h  without 
geometrical  effects 

Analytical 
h  =  O 

0.1 

4.232  x  10  ' 

4.21  1  x  10  ‘ 

3.800x10  ' 

3.741  x  l(t  ' 

2.346  x  |0  ' 

2.4 1 4  x  10  ' 

2.671  x  10  < 

0.5 

2.212  x  10  ; 

2.201  x  10  ; 

2. 1 42  x  1 0  : 

2.067  x  10  ; 

1. 1  264  x  ID  ; 

1.31  3  x  10  ■ 

1.432x10  - 

1.0 

0.2339 

0.2328 

0.2443 

0.2355 

0.1408 

0.1482 

0.1585 

1.5 

0.7504 

0.7469 

0.7559 

0.7234 

0.4675 

0.4874 

0.5256 

2.0 

1.5574 

1.5474 

1.41 13 

1.3379 

0.9831 

1.0029 

1.1  10 
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Table  1  shows  the  comparison  between  the 
results  of  various  approximations  and  the  exact 
calculations  for  a  spherical  inclusion.  The  prop¬ 
erties  of  the  layer  have  been  taken  to  be  the 
mean  of  those  of  the  inclusion  and  matrix.  The 
thickness  h  is  taken  to  be  0.1«,  a  being  the  radius 
of  the  sphere.  It  is  seen  from  this  table  that  both 
geometrical  and  inertial  effects  are  important  in 
the  determination  of  the  scattering  cross-section. 
Other  conclusions  that  can  be  drawn  from  this 
table  are  as  follows. 

( 1 )  Full  order  h  predictions  are  pretty  close  to 
the  exact  results  for  low  and  intermediate 
frequencies. 

(2)  A  simplified  approximation  obtained  by 
keeping  only  the  inertial  effects  is  fairly  good. 

(3)  The  approximations  resulting  from  keep¬ 
ing  the  jump  in  u  (as  is  done  commonly)  is  inferior 
to  both  the  membrane  and  the  inertial  approxi¬ 
mations.  In  fact,  this  results  in  a  prediction  that  is 
lower  than  that  for  an  inclusion  without  an  inter¬ 
face  layer.  It  is  found  that  the  prediction  of  an 
O  (h)  aproximation  without  the  geometrical 
effects  is  also  lower. 

(4)  The  presence  of  the  interface  layer  around 
the  particle  changes  the  scattering  cross-section 
measurably.  The  implication  of  this  for  a  particle- 
reinforced  composite  material  would  be  to 
change  the  phase  velocity  and  attenuation  of 
ultrasonic  waves.  Thus,  measurements  of  these 
quantities  would  allow  characterization  of  the 
interface.  However,  our  study  suggests  that,  in 
order  to  carry  out  the  inversion  correctly,  it  is 
necessary  to  include  the  effects  of  curvature  and 
inertia  of  the  layer.  A  simple  spring  model 
neglecting  these  effects  may  lead  to  erroneous 
conclusions. 

In  order  to  demonstrate  the  inertial  effects  we 
present  in  Fig.  1  a  comparison  between  the  pre¬ 
dicted  scattering  cross-sections  of  a  prolate 
spheroid  with  an  aspect  ratio  bja  of  2  to  1  for  a 
longitudinal  wave  incident  at  45°  to  the  symmetry 
axis.  These  tesults  are  obtained  using  the  mem¬ 
brane-type  approximation  and  for  the  following 
three  types  of  interface  layers:  type  1  is 


_  A  +  A, 

"  “  2 

P+P, 

Al  =  ~ 

P+P  1 

P"~  2 

type  2  is 

An  +  2/r0=  36.8  GPa 
A>  =  8.9  GPa 

P»  =  P 
and  type  3  is 
A(l  +  2ju„  =  36.8  GPa 
=  8.9GPa 
p„  =  1-353  kgm 

It  is  seen  from  these  values  that  as  the  density 
decreases,  so  does  the  scattering  cross-section. 
The  substantial  increase  in  the  scattering  cross- 
section  in  the  presence  of  the  interface  layer  is 
also  noticeable. 

4.2.  Effect  of  an  interface  layer  on  dispersion  of 
guided  waves  in  a  laminated  plate 

We  next  present  some  results  showing  the 
measurable  changes  in  the  phase  velocity  disper- 


0.0  0.5  1,0  Kb  1.5  2.0 

Fig.  1 .  Total  scattering  cross-section  of  a  prolate  f  bja  =  2/ 1 ) 
spheroid  for  an  incident  P  wave  propagating  at  45°  with  the 
symmetry  axis. 


T  ABLE  2  Properties  of  0°  and  90°  laminae  and  of  the  interface  layer 


<x'l()"Nm  (  xlO"  N  m  ;)  (x‘l()"Nm  -)  (x)o"Nm  -’)  jx'i()"Nm 


0°  lamina  1.0673  0.1392  0.0644  0.0350  0.0707 

interface  0.0865  0.0865  0.0475  0.0195  0.0!  95 

90°  lamina  0.1392  1.0673  0.0644  0.0707  0.0350 
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sion  in  a  cross-ply  laminated  plate  caused  by 
interface  bond  layers.  We  consider  graphite-fiber- 
reinforced  laminae.  The  effective  elastic  proper¬ 
ties  of  a  lamina  and  the  interface  bond  layers  are 
listed  in  Table  2. 

Figures  2-7  show  the  variations  in  phase 
velocities  with  frequency  in  a  laminated  plate  with 
19  laminae  with  and  without  interface  layers. 
Figures  2  and  3  are  for  propagations  along  x  (0°) 
and  y  (90°)  directions  respectively  when  there  are 
no  interface  layers.  Figures  4  and  5  show  the 
effect  of  interface  layers.  The  interface  layer 
thickness  is  taken  to  be  one  tenth  of  a  lamina 
thickness  and  the  ratio  of  the  densities  of  a  layer 
and  a  lamina  is  0.67.  To  show  the  effect  of  density 
(inertia)  we  show  in  Figs.  6  and  7  the  dispersion 


NORMALIZED  EREOuENCt 

Fig.  2.  Dispersion  of  guided  waves  in  a  ()°/9()°/()'... 
laminated  Id  layers)  plate  for  propagation  in  the  0°  direc¬ 
tion.  - ,  symmetric  modes; - ,  antisymmetric  modes;  no 

interface. 


Fig.  3.  Same  as  in  Fig.  2  for  propagation  in  the  90°  direction. 


curves  when  the  density  ratio  is  1.5.  In  all  these 
figures  the  non-dimensional  phase  velocity  and 
frequency  are  defined  as 

tu/k 

HkwpwF 

Q  -  (°H 

~  2K  cwpfcTF 

where  H  is  the  total  thickness  of  the  plate. 

The  following  conclusions  can  be  drawn  from 
these  figures: 

( 1 )  It  is  seen  from  Figs.  2  and  3  that  at  low  fre¬ 
quencies  the  dispersion  for  propagation  along  0° 

1  5 - - - 

1  3  : 
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Fig.  4.  Dispersion  of  guided  waves  in  a  079070°... 
laminated  plate  ( 19  laminae)  with  interface  layers  between  0° 
and  90°  laminae.  The  ratios  of  the  thicknesses  and  densities 
of  a  lamina  and  a  layer  are  10  and  1.5  respectively.  Propaga¬ 
tion  is  in  the  0°  direction. - .  symmetric  modes; - .  anti¬ 

symmetric  modes. 
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Fig.  5.  Same  as  in  Fig.  4  for  propagation  in  the  90°  direction. 
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Fig.  6.  The  configuration  is  the  same  as  in  Fig.  4,  hut  the 
density  ratio  is  0.67. 
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Klu.  7.  Same  as  in  l  ip.  6  for  propagation  in  the  90°  direction. 


and  40°  directions  is  the  same.  Thus  the  plate 
behaves  quasi-isottopically  at  low  frequencies. 
However,  at  hitjh  frequencies  the  dispersions  for 
propagation  in  these  two  directions  are  quite 
different. 

(2)  Comparison  of  Figs.  4  and  5  with  Figs.  2 
and  3  shows  that  the  presence  of  the  interface 
layers  lowers  the  phase  velocities  and  the  cut-off 
frequencies.  The  effect  is  amplified  at  high 
frequencies. 

(’3)  Finally,  Figs.  6  and  7  show  that  raising  the 
density  of  the  interface  layers  does  not  change  the 
dispersion  of  the  first  few  modes  appreciably. 
However,  with  increasing  frequency  the  differ¬ 
ences  between  F  igs.  4  and  5  and  Figs.  6  and  7 
become  appreciable.  At  certain  frequency  bands 


these  differences  are  quite  large.  This  suggests 
that,  in  order  to  extract  the  bond  layer  properties, 
judicious  choice  of  frequency  and  modes  can  be 
made  to  obtain  optimum  results. 

5.  Conclusion 

The  effect  of  interface  and  interphase  layers  on 
wave  propagation  in  a  composite  medium  has 
been  studied  in  this  paper.  Two  problems  have 
been  considered:  scattering  by  a  reinforcement 
particle  and  guided  waves  in  a  laminated  plate. 
The  results  of  this  investigation  show  that  in  a 
particulate-reinforced  composite  the  presence  of 
an  interface  layer  between  the  particles  and  the 
matrix  would  cause  significant  changes  in  the 
damping  of  ultrasonic  waves.  Also,  phase  velocity 
dispersion  in  a  laminated  plate  changes 
measurably  in  the  presence  of  interface  bond 
layers.  It  is  found  in  both  cases  that  the  density  of 
the  layer  plays  an  important  role  in  modifying  the 
wave  propagation  characteristics.  In  addition,  for 
curved  interfaces  the  curvature  effect  should  be 
taken  into  account.  Thus,  in  order  to  characterize 
interfaces  properly,  it  would  be  necessary  to  take 
into  a<  count  properly  not  only  the  stiffnesses  but 
also  the  inertial  and  curvature  effects. 
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